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Abstract 

A one-loop correction of the quasilocal energy in the Schwarzschild back- 
ground, with flat space as a reference metric, is performed by means of a 
variational procedure in the Hamiltonian framework. We examine the gravi- 
ton sector in momentum space, in the lowest possible state. An application to 
the black hole pair creation via the Casimir energy is presented. Implications 
on the foam-like scenario are discussed. 



I. INTRODUCTION 

It is generally accepted that if Einstein gravity has a ground state, this is represented 
by Flat space. Many attempts to discover a possible lower minimum different from the Flat 
one have been carried out and it seems nowadays that, at least at classical level, the lowest 
energy configuration state with E = is attributed to Flat space IQ. From the quantum 
mechanical point of view we can obtain information on the stability of the spacetime by 



1 



means of saddle point methods: if a lower state exists, this can be reached by quantum 
tunneling. As a consequence a negative mode will appear in the Lichnerowicz operator. 
Actually a tunneling process was discovered by Gross, Perry and Yaffe 0], but it involves 
a flat space having a temperature T 7^ 0, which means that we have not considered as 
initial state the correct "vacMum" . Thus also this result seems to corroborate the stability 
of fiat space under quantum fluctuation via nucleation of black holes. Recently a different 
mechanism involving both a topology change and a black hole pair production has been 
considered Q. In this paper we are investigating the possibility of computing such effects 
in a Hamiltonian approach in presence of quasilocal energy. To this purpose we begin by 
fixing the background manifold M. which is represented by the Schwarzschild spacetime. In 
particular we are interested in a constant time section of M., best known as Einstein- Rosen 
bridge with wormhole topology x S^, defining a bifurcation surface which divides S into 
two parts denoted by S+ and S_. Our purpose is to consider perturbations at S, which 
naturally define quantum fluctuations of the Einstein- Rosen bridge and evaluate quasilocal 
energy corrections to one-loop with the assumption of neglecting quantum excitations on the 
boundaries, motivated by the fact that in asymptotic spacehke directions, quasilocal energy 
approaches the AVM. term defined in the asymptotic region, whose quantum fluctuations 
are unphysical. Nevertheless this approximation is not valid when we consider finite distance 
located boundaries and multi-wormholes configuration 0,3. In this context, we will examine 
the effect of perturbations on the possible ^''ground state" by means of a variational approach 
applied on gaussian wave functional. Indeed, if we discover the existence of negative modes 
in this approximation, it is quite reasonable to think of a tunneling process which moves 
spacetime from the false vacuum to the true one. The rest of the paper is structured as 
follows, in section we have borrowed from Refs. p,^] the general expressions for Hamil- 
tonian and quasilocal energy, in section |ITT|, we analyze the stationary Schrodinger equation 
coming from the perturbed wormhole metric and we give some of the basic rules to perform 
the functional integration for the Hamiltonian approximated to second order, in section 
we analyze the spin-2 operator or the operator acting on transverse traceless tensors. We 
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summarize and conclude in section 0. 



II. QUASILOCAL ENERGY AND ENERGY DENSITY CALCULATION IN 
SCHRODINGER REPRESENTATION 

Although it is not necessary for the forthcoming discussions, let us consider the maximal 
analytic extension of the Schwarzschild metric, i.e., the Kruskal manifold whose spatial slices 
E represent Einstein- Rosen bridges with wormhole topology S*^ x R^. Following Ref. [§, 
the complete manifold Ad can be taken as a model for an eternal black hole composed of 
two wedges Ai^ and located in the right and left sectors of a Kruskal diagram. The 
hypersurface S is divided in two parts S+ and S_ by a bifurcation two-surface Sq. On E 
we can write the gravitational Hamiltonian 

Hp = H -Ho = d^xiNH+N'Hi) 
+ - [ d\N^ (k-k^)- - [ d\N^ (k - k^) , (1) 

K Js+ ^ ^ K JS- ^ ^ 

where k = SttG. The Hamiltonian has both volume and boundary contributions. The 
volume part involves the Hamiltonian and momentum constraints 

n = (2k) G,,ki'^^^ - J^gR/ {2k) = 0, 

= -2vr5^. = 0, (2) 

where Gijki = \ {gikdji + diWjk — gijQki) is the supermetric and R denotes the scalar curvature 
of the surface S. The volume part of the Hamiltonian (0) is zero when the Hamiltonian 
and momentum constraints are imposed. However, for the fiat and the Schwarzschild space, 
constraints are immediately satisfied, then in this context the total Hamiltonian reduces to 

Hp = - f d^xN^ [k-k°)- - f d^xN^ {k - k°) . (3) 
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Quasilocal energy is defined as the value of the Hamiltonian that generates unit time trans- 
lations orthogonal to the two-dimensional boundaries, i.e. 



Etot — -f^quasilocal — — — — E_, 

E+^-f d^x^/d (k - k^) 

K JS+ ^ ' 



K 



f (fx^^k-k'^). (4) 



where |A^| = 1 at both 5*+ and S_. E^ot is the quasilocal energy of a spacelike hypersurface 
E = E+ U E_ bounded by two boundaries ^5"+ and ^S- located in the two disconnected 
regions M+ and M_ respectively. We have included the subtraction terms for the energy. 

represents the trace of the extrinsic curvature corresponding to embedding in the two- 
dimensional boundaries ^5"+ and ^S- in three-dimensional Euchdean space. Following the 
same scheme of the boundary subtraction procedure, we would like to discuss the possibility 
of generalizing such a procedure. To this end, by looking at the Hamiltonian structure, we 
see that there are two classical constraints 

= n 

(5) 

which are satisfied both by the Schwarzschild and Flat metric and two quantum constraints 

(6) 

W"^ = 

= is known as the Wheeler- De Witt equa,tion (WDW). Nevertheless, we are interested 
in assigning a meaning to 




^Schw. ^Flat 



|gqua^ilocal| ^) 



where \E' is a wave functional whose structure will be determined later and ifSchw. |^f/'Fiat^ 
is the total Hamiltonian referred to the different spacetimes. Note that the first term of 



(0) is simply the extension of the quasilocal energy subtraction procedure generahzed to the 
volume term. Note also that if the expectation value is calculated on the wave functional 
solution of the WDW equation, we obtain only the boundary contribution. Then to give 
meaning to (0), we adopt the semiclassical strategy of the WKB expansion. By observing 
that the kinetic part of the Super Hamiltonian is quadratic in the momenta, we expand the 
three-scalar curvature / cPx^^gR^^'' up to quadratic order and we get 



d X 



(8) 



where h is the trace of hij. On the other hand, following the usual WKB expansion, we 
will consider \1/ ^ Cexp(z5'). In this context, the approximated wave functional will be 
substituted by a trial wave functional according to the variational approach we would like 
to implement as regards this problem. 

III. GAUSSIAN WAVE FUNCTIONAL AND ENERGY DENSITY 
CALCULATION IN SCHRODINGER REPRESENTATION 

To actually make such calculations, we need an orthogonal decomposition for both and 
to disentangle gauge modes from physical deformations. We define the inner product 

(/i, k) := ^G'^^'h,^ (x) hi (x) d'x, (9) 

by means of the inverse WDW metric Gijki, to have a metric on the space of deformations, 
i.e. a quadratic form on the tangent space at h, with 

Qijkl ^ (^gikgjl gilgjk _ 2^*^^'='). (lO) 

The inverse metric is defined on co-tangent space and it assumes the form 

(P, q) ■■= I VaG^jkip'' (x) q'' {x) d'x, (11) 

so that 
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(12) 



Note that in this scheme the "inverse metric" is actually the WDW metric defined on 
phase space. Now, we have the desired decomposition on the tangent space of 3-metric 
deformations 



hij = -hg,j + {LC)^. + hf^ 



where the operator L maps into symmetric tracefree tensors 



(13) 



(LO., = v.O + - ^g^, (V ■ 



Then the inner product between three-geometries becomes 



(/i, h) := / ^G'^^^hij (x) hki (x) £x 

J X] 



(14) 
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With the orthogonal decomposition in hand we can define a " Vacuum Trial State" 



^^J[h,^{lf)]=^fexp\-— 



x,y 



hK-^h 



Trace 



x,y 



(15) 



, (16) 



which will be used as a probe for the gravitational ground state. This particular expression 
is useful because the functional can be represented as a product of three functional defined 
on the decomposed tensor field 



vj/ [h,, (x^)] = m ht^ (1^) ^ (LO,, ^ -9^Jh {Ig 
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(17) 



h:^j is the tracefree-transverse part of the 3D quantum field, {L^)ij is the longitudinal part 
and finally h is the trace part of the same field. {■,-)xy denotes space integration and 
K"^ is the inverse propagator containing variational parameters. The main reason for a 
similar "^nsat^" comes from the observation that the quadratic part in the momenta of the 



Hamiltonian decouples in the same way of eq.(|T5|). Note that the decomposition related to 



the momenta is independent of the choice of the functionaL To calculate the energy density, 
we need to know the action of some basic operators on [hij] . The action of the operator 
hij on |\[^) = \E' [hij] is realized by 

hij{x)\^) = hij{lt)^[hij]. (18) 

The action of the operator TCij on |\E'), in general, is 

Tr,^(x)\^) = -tjj^-^^[h,^]. (19) 

The inner product is defined by the functional integration: 

(^1 I ^2) = / [Vhi,] {h,,} ^2 {hki} , (20) 

and the energy eigenstates satisfy the stationary Schrodinger equation: 

J d^xH i^-tjj^^^, ( a^)| ^ {K,} = E^> {K,} , (21) 

where ^{^—i ^^ ^(x) ' (••^) } Hamiltonian density. Note that the previous equation 

in the general context of Einstein gravity is devoid of meaning, because of the constraints. 
However in the semiclassical context, we can give a meaning to eq. (PTf) , where a semiclassical 
time is introduced in the same manner of Refs. PJT(]||. There, a Schrodinger equation of the 
form 



I- 



i7|2^^ (22) 



dt 

is recovered by the WDW equation approximated to second order for a perturbed minisu- 
perspace Friedmann model without boundary terms. When asymptotically flat boundary 
terms are present we have to take account of such contributions in the WKB expansion 
such as in Ref. . However in this paper only gravitational transverse-traceless modes are 
considered on the fixed curved background and is substituted by a trial wave functional. 
To further proceed, instead of solving (PH), which is of course impossible, we can formulate 
the same problem by means of a variational principle. We demand that 



/ 






/ 







(23) 



be stationary against arbitrary variations of ^{/ij^}. The form of \ H\2 \ '^j can be 
computed as follows. We define normalized mean values 



/ 




S<Pxgi^{lt)\^>[gi^} P 


/ 







(24) 



ah (^) ail it) + k^.^i t) 



(25) 



/ 


Vgi,_ 


I d^xgij ( 


-^)ahm\'^[aii} ? 


/ 


-Dgi, 


\^{at^}? 



It follows that, by defining h^, = gij — gij, we have 



Vh 



hij{lf)\^{hi-^+gl-^}\'=0 



(26) 



(27) 



and 



Vh 



d^xhi^ (If) 4 it) I ^ {hij + g^j} 



Vh 



"^H + ai} 



(28) 



Nevertheless, the application of the variational principal on arbitrary wave functional does 
not improve the situation described by the eg. (pT]) . To this purpose, we give to the trial 
wave functional the form 



h 



^^-p[-^{ia-a)K-Ha-a)]^^^ 



(29) 



We immediately conclude that 



^|7r,:^(^) 1^) = 



(30) 



where tt^^ is the TT momentum. In Appendix 0, we will show that 
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(^) 4 it) 1^) = \Kr^M V) • (31) 

Choice (^) is related to the form of the Hamiltonian approximated to quadratic order in 
the metric deformations. Indeed, up to this order we have a harmonic oscillator whose 
ground state has a Gaussian form. By means of decomposition (|T3p, we extract the TT 
sector contribution in the previous expression. Moreover, the functional representation (|T^ 
eliminates every interaction between gauge and the other terms. Then for the TT sector 
(spin-two), one gets 

d'x^R'^'^ ^^J^ d'x^ [h^'^ (A2); hi - 2hR,^h^'=] , (32) 

where (A2)" := + 2i?". The latter term disappears because the gaussian integration 

does not mix the components. Then by collecting together eq.(^) and eg . ([3l|) , one obtains 
the one-loop-like Hamiltonian form for TT deformations 

H^ = ^i d'x^G^^'^ [k-^^ (x, x\^,, + ( A2); (x, x),,,J . (33) 

The propagator (x, x)^^^; comes from a functional integration and it can be represented 

as 

where ( af) are the eigenfunctions of and Xn ip) are infinite variational parameters. 



IV. THE SPECTRUM OF THE SPIN-2 OPERATOR AND THE EVALUATION 
OF THE ENERGY DENSITY IN MOMENTUM SPACE 

The Spin-two operator is defined by 

(A2); := -A6; + 2R] (35) 

where A is the curved Laplacian (Laplace-Beltrami operator) on a Schwarzschild background 
and Rj is the mixed Ricci tensor whose components are: 



f —2m m m 
= diag { — — , — 



(36) 



where 2m = 2MG. This operator is similar to the Lichnerowicz operator provided that we 
substitute the Riemann tensor with the Ricci tensor. This is essentially due to the fact that 
the Riemann tensor in three-dimensions is a linear combination of the Ricci tensor. In (|37D 
the Ricci tensor acts as a potential on the space of TT tensors; for this reason we are led to 
study the following eigenvalue equation 



[-A6] + 2R';) h\ = E^h] 



(37) 



where E'^ is the eigenvalue of the corresponding equation. In doing so, we follow Regge 
and Wheeler in analyzing the equation as modes of definite frequency, angular momentum 
and parity. We can specialize to the case with M = without altering the contribution 
to the total energy because of the spherical symmetry of the problem. We recall that L 
is the quantum number corresponding to the square of angular momentum and M is the 
quantum number corresponding to the projection of the angular momentum on the z-axis. 
In this case, Regge- Wheeler decomposition [n] shows that the even-parity three-dimensional 
perturbation is 



diag 



9m \ 



Yio (^,0) 



(38) 



Representation (^) shows a gravitational perturbation decoupling. For a generic value of 
the angular momentum L, one gets 

-AiH{r)-^H{r) = EfH{r) 



-AiK{r) + ^K{r) = EfK (r) 

-AiK{r) + '^K{r)=EfK (r) . 
The Laplacian restricted to S can be written as 

2m\ f2r-3m\ d / (/ + 1) 



(39) 



A 



r / dr^ 



dr 



(40) 
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Defining reduced fields 



H(r) 



h{r) 



Kir] 



k (r) 



(41) 



and passing to the proper geodesic distance from the throat of the bridge defined by 

dr 



dx 



1 - 



2M 



whose integrated form is 



x = 2m{\ J l + ln W hW 1 

2mV2m VV 2m V 2m 



the system (|3^) becomes 



-£,h {%) - (x) h (x) = Efh (x) 



< -£,k{x) + V+ {x)k{x) = Efk{x) 



k (x) + V+ (x) k (x) = Efk ix) 



with 



V ix) = \ , / =F — 



r^fx) ' r(xf 



(42) 



(43) 



(44) 



(45) 



This new variable represents the proper geodesic distance from the wormhole throat such 
that 

when r — > oo, x ~ r , (x) — > 



when r — > tq, x ~ 0, (x) 



/(/ + 1) 3m 

5 — =F = const, 



(46) 



'0 '0 

where tq satisfies the condition tq > 2m. The solution of (^if), in both cases (fiat and curved 
one) is the spherical Bessel function of the first kind 



Jo (px) = \/-sin (px) 



(47) 
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This choice is dictated by the requirement that 



h{x) ,k (x) when x — > (alternatively r — > 2m) . 



(48) 



Then 



X _ Jo (px) jo (py) J_ 



(49) 

Substituting (|49| ) in (|33| ) one gets (after normalization in spin space and after a rescaling of 
the fields in such a way as to absorb 1^) 



where 



^1,2 (P> - P + -2 T -3-, 



(50) 



(51) 

'0 '0 

Aj (p) are variational parameters corresponding to the eigenvalues for a (graviton) spin-two 

particle in an external field and V is the volume of the system. 

1 

By minimizing ( ^OD with respect to Aj (p) one obtains Aj (p) = [Ef {p, m, /)] ^ and 

'(m,A) = — EE/ dp2^Ef{p,m) 
;=0 i=l 



E 



(52) 



with 



with p'^ + 



/(/ + !) 3m 



Thus, in presence of the curved background, we get 



Eim) = —-J2j^ dpp'[^Jp^ + cl + ^p^ + c\^ 



(53) 



where 



2 _ / (/ + 1) 3m 

~ ^^2 T —3-' 



while when we refer to the flat space, in the spirit of the subtraction procedure, we have 



m = and = Then 

^'0 
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V 1 °° r°° / I \ 



(54) 



Now, we are in position to compute the difference between ( ^3]) and (^41). Since we are 
interested in the UV hmit, we have 



AE (m) = E{m)- E (0) 



y \ ^ i-oo 



+ c_ + Jp^ + ci — 2wp^ + c 



Y I °° roo 



and for >> c^,c^, we obtain 



V 1 
2^2 



E / dpp^ 



+ 



2 (p. 



sip) ^^^2 (p ) sip) 



(55) 



-2 - 



4 L 



y /■°° dp 



27r2 8 Jo J9 



We will use a cut-off A to keep under control the UV divergence 



dp 
p 



c doc 



X 



In 



A 



where A < mp. Thus AE (m) for high momenta becomes 



AE (m) 



V f3m\ 1 j^/^oA^^ 
16 I 3m 



(56) 



(57) 



(58) 



27r2 16 

Remark It is known that at one-loop level Gravity is renormalizable only in flat space. In a 
dimensional regularization scheme its contribution to the action is, on shell, proportional to 
the Euler character of the manifold that is nonzero for the Schwarzschild instanton. Although 
in our approach we are working with sections of the original manifold to deal with these 
divergences one must introduce a regulator that indeed appears in the contribution of energy 



density, which we presume should be taken of order the Planck mass mp ~ V in the 
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same spirit of Ref. 0]. This agrees with Wheeler's point of view concerning the gravitational 



vacuum fluctuations, where a natural cut-off is introduced: the Planck length /p ~ v G [|T5 
At this point we can compute the total energy, namely the classical contribution plus the 
quantum correction up to second order. Recalling the expression for quasilocal energy 

^quasilocal = " / d^Xy^ (k - k^) - - (fx^ [k - , (59) 

and by using the expression of the trace 



k 



Vh 

we obtain at either boundary that 

-2r 



L(VA„.) (60) 



k = (61) 



r 



where we have assumed that the function r,y is positive for and negative for S^. The 
trace associated with the subtraction term is taken to be = — 2/r for and k^ = 2/r 
for B_. Then the quasilocal energy with subtraction terms included is 

^quasilocal = - = {v [1 - |r„ |])^^^^ - (r [1 - |r„ |])^^^_ . (62) 

Note that the total quasilocal energy is zero for boundary conditions symmetric with respect 



to the bifurcation surface 5*0 . The energy E^ (E-) tends to the AVM. mass M [|12| whenever 
the boundary {BJ) tends to right (left) spatial infinity. As an illustration, consider the 
case when the boundary B^ is located at right-hand infinity (?/+ = +00) and the boundary 
B^ is located at The total energy for the stable modes is 



M -r 



V fSMGy 1 , /r^A^ \ 



2n 



while the total square energy for the unstable modes is 

.2 



E 



2 



a 



{MGf 



with = 0, 242 and will be computed in the appendix |C[ Then the one loop total energy is 
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M -r 



2MG\ 



-) 



V f3MG\ 1 , /rM^ 
I 1 In ' 



(64) 



27r2 \^ J 16'"\3MGJ ' MG' 
If we consider only one of the wedges, we obtain the correction relative to the positive 
(negative) AT>Ai mass M , depending on the location of the wedge. One can observe that 



AE (M) oo when M ^ 0, for tq = 2GM 



(65) 



and 



AE (M) ^ when M ^ 0, for tq 7^ 2GM. 



(66) 



Note that this singular behaviour is independent of boundary conditions since it is related 
to the volume term. 



V. SUMMARY AND CONCLUSIONS 

We started from the problem of defining (semiclassical) quantum corrections to the 
quasilocal energy. By means of a variational approach with Gaussian wave functionals, 
an attempt to calculate such a correction was made. By construction, we used the subtrac- 
tion procedure given in Refs. to avoid divergences coming from boundaries. Despite 
the constraint equations, this calculation is based on an extension of the subtraction proce- 
dure involving volume terms in the semiclassical regime. Excitations coming from boundary 
terms have been neglected to avoid the unphysical situation of having contributions deriving 
from infinity. In this context, the extended subtraction procedure corresponds to the differ- 
ence between zero point energies calculated in an asymptotically flat background referring 
to a fiat background. This procedure eliminates the UV divergence of the free gravitons, 
leaving the contribution of the curved background related to an imposed by hand UV cut-off. 
Note that the subtraction procedure at one loop has the correct ingredients to be related 
to the Casimir energy. This seems to give some information about the vacuum behaviour. 
Indeed, if we look at symmetric boundary conditions with respect to the bifurcation surface 



5*0, then eq. (|64D can be interpreted as an energy gap measuring the probability of creating a 
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black hole pair by a topological fluctuation; it also indicates that flat spacetime is unstable 
with respect to pair creation. In the introduction we mentioned the impossibility of gen- 
erating single black holes by quantum fluctuation in a flat space having zero temperature, 
because the energy would not be conserved [0,|16|. However, in the neutral black hole pair 
scenario [Q, where each component resides in a different universe, the energy is conserved 
provided that the boundaries are symmetric with respect to the bifurcation Sq. Note that 
this calculation can be related to the quantity 

r -Pi— hole -pBlackHolePair t r'<-r\ 

1-hole = ^ ^ . (67) 

-Tflat -fflat 

To generalize a little more, suppose to enlarge this process from one pair to a large but flxed 
number of such pairs, say n. What we obtain is a multiply connected spacetime with n holes 
inside the manifold, each of them acting as a single bifurcation surface with the sole condition 
of having symmetry with respect to the bifurcation surface even at flnite distance. Let us 
suppose the interaction between the holes can be neglected, i.e., let us suppose that the 
total energy contribution is realized with a coherent summation process. This is equivalent 
to saying that the wave functional support (here, the semiclassical WDW functional) has a 
finite size depending only on the number of the holes inside the spacetime. It is clear that 
the number of such holes cannot be arbitrary, but is to be related with a minimum size of 
Planck's order. Thus, assuming a coherency property of the wave functional and therefore 
of the A^-holes spacetime, eq.(B^ has to be generalized to 



m -Pn— holes -Pfoam ir'n\ 

TN-hoies = —E, s — • (68) 

-iflat -fflat 



Recall that a hole, here, has to be understood as a wormhole. Then, what eq.(|68|) suggests 



is the possibility of generating a foamy spacetime, with wormholes as building blocks. 



APPENDIX A: CONVENTIONS AND SCALAR CURVATURE EXPANSION 



Riemann tensor 
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ijm Tni,j ji,m ' ja mi ma ji 



(Al) 



Ricci tensor 



D pi 



(A2) 



• Scalar curvature 



(A3) 



• In three dimensions, the Weyl tensor vanishes, then the Riemann tensor is completely 
determined by Ricci tensor 



Rlijm — 9ljRim QlmRij Qij^lm + 9imRlj- 



(A4) 



• Second order scalar curvature 



/ 



(fx 



2' 



1 

-I 
2 



(A5) 



APPENDIX B: THE KINETIC TERM 

The Schrodinger picture representation of the kinetic term is 

i - ki ( \ 



dhij (x) dhki {x) J ' 

We have to apply this quantity to the gaussian wave functional \^). This means that 

6^^ [h] 



(Bl) 



tt'^ (x) tt^^' (x) I*) 



6hij (x) 6hki {x) 
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■K-'^^^^^'"'"\xy)h,,n"{y")'^[h]. (B2) 

By functional integration 

\hkn' {y') hk"i" {y")\ ^) = K^k'i'wn iy',y") (^1^) • (B3) 



Then 



becomes 



^ n'^ (x) Tx^^ ix) ^ 



1^-1 Wfo) (3,^^) 



K(k'i')(k"i") {y',y") (^1^) 

= ^K-'^mj) (^,^) (^y^)' (vl/|vl/) . (B4) 
Then the expectation value of the kinetic term, with the Planck length reinserted, is 

(T) = ^ I d'x^ (G.^mK-'^"'^^'^^ (x, x)) , (B5) 

APPENDIX C: SEARCHING FOR NEGATIVE MODES 

In this paragraph we look for negative modes of the eigenvalue equation (|37D. For this 
purpose we restrict the analysis to the S wave. Indeed, in this state the centrifugal term 
is absent and this gives the function V {x) a potential well form, which is different when 
/ > 1, where / is the angular momentum. However, this potential form is valid only for the 
H component, i.e. 

-AH (r) - (r) = -E^H (r) , 



with E"^ > 0. Passing to the reduced field H (r) = we obtain 

d I / 2mdh\ /—3m h 
'l ^ I + (— — + 



dr \ \ r dr I \ J ^\ 2m 

Making the substitution (|43D , the equation becomes 



dx d I I 2m dh dx 
dr dx W r dx dr 



d / dh\ / 3m ^2 \ ^ g 
dx \ dx \ r^ 



where ^ 

dr 



2m 



— ^==. Near the horizon x ~ 4my^2^ — lor r ^ 2m ^1 + j • Define 

, 2 



p = 1 + (^ij with y = Then we get 



d fdh\ 
dy \dy) 



3m 



(2m) p3(?/) 



d'^h 
dy' 



+ 



V 



2(1+1 



+ A 



/i = 0, 



where A = (2m) E"^. Expanding the potential around y = 0, one gets 

rf^/i / 3 / 3 



dy^ 



(i - -r 

2 V 4:^ 



\]h = 0. 



We see that this is a quantum harmonic oscillator equation whose spectrum is well known, 
that is En = huj (^n+ . However, we have to remark that the original spectrum must be 
contained in the interval , where — | represents the bottom of the approximated potential. 



Then En = l + K 



I m + I j , where u) 



and h = 1 in natural units. Thus 



3^2 
2^ 



n 



We see that 



Ao = -0.975 and Ai = — (-0.06) 

2 



Ai 



There is only one eigenvalue. The same problem can be approached with the Rayleigh-Ritz 
method along a numerical integration and the result is Aq = —1.094. 
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